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Abstract 



I examine the existence of zero energy fermion solutions (zero modes) on cosmic 
strings in an SO{10) grand unified theory. The current carrying capabihty of a cosmic 
string formed at one phase transition can be modified at subsequent phase transitions. 
I show that the zero modes may be destroyed and the conductivity of the string altered. 
I discuss the cosmological implications of this, and show that it allows vorton bounds 
' to be relaxed. 

^ ' 1 Introduction 

Q^l Cosmic strings, which are a type of topological defect, arise in many grand unified theories. 



Large quantities of them may be produced at phase transitions in the early universe. A net- 
work of cosmic strings could explain the observed anisotropy in the microwave background 
radiation and the large scale structure of the universe |||]. 
^ ' It has been realised in the past few years that cosmic strings have a far richer microstruc- 

. ture then previously thought In particular, the presence of conserved currents in the 

spectrum of a cosmic string has profound implications for the cosmology of the defects. In 
this paper I examine currents with fermion charge carriers, particularly massless ones. 

Currents may provide a method of detecting strings. If they form at high energy scales, 
and are charged, the resulting electromagnetic field may be detectable [^. Decaying currents 
may explain observed high-energy cosmic rays ||3|, Q], and could also provide a mechanism 
for baryogenesis js], There are several processes which can create currents on strings. 
These include interaction with the plasma and collisions between cosmic strings. Charged 
currents can also be generated by magnetic or electric fields. Stable relics, vortons, can 
form as collapsing string loops are stabilized by the angular momentum of the trapped 
charge carriers Q. If they form at high energy scales, these relics can dramatically alter 
the evolution of string networks. Vortons formed at low energies may provide a dark matter 
candidate [p. 

As the universe expands, a network of cosmic strings will stretch, suggesting that its en- 
ergy density will grow relative to everything else. This would lead to the universe becoming 
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string-dominated, which is strongly ruled out by observation. More careful analysis shows 
that loops of string will break off from the network. These then contract, losing energy 
by gravitational radiation. This avoids string domination. If the loops are stabilized by 
conserved currents, this mechanism fails, allowing the corresponding theory to be severely 
constrained 

The cosmological implications of the vortons are most pronounced when the universe 
has become matter-dominated. If they decay during the era of radiation domination, the 
cosmological catastrophe may be avoided and the vorton bounds evaded. It has recently 
been realized that subsequent phase transitions can have a considerable effect on the micro- 
physics of cosmic strings. Massless fermion currents on cosmic strings can be both created 



and destroyed [10|. In this paper I examine the destruction of such currents in grand 
unified theories with an 50(10) symmetry group. 

I describe a simple cosmic string model in Section and show that it has zero energy 
fermion solutions, called zero modes. This result was originally derived by Jackiw and 
Rossi I then show how the zero-mode solutions can be extended to massless currents. 

In Sectional discuss the fate of zero modes on strings formed in an SO{10) grand unified 
symmetry At high temperatures the theory resembles the toy model of Section 0. I 
investigate the implications of the electroweak phase transition for the zero modes, and 
show that they do not survive it. 

I consider the implications of spectral flow in Section ^, and deduce that the currents 
acquire a small mass. The string current can now dissipate. This allows vortons to decay 
and weakens the cosmological bounds on such models [^. An important feature that allows 
zero modes to be removed is the presence of a massless particle that mixes with the zero 
mode after the transition. Finally, I summarize the conclusions. 

2 Fermion Zero Modes on an Abelian Cosmic String 

The Abelian Higgs model has the Lagrangian 

£ = (D.^riD^cP) - \f,,F^" - - rj'f (1) 

with F^i, = dfj,Ay — dyAfj_. As well as the usual (p = const, solution, it has cosmic string 
solutions of the form 

^ = rif{r)e'-o (2) 

A, = n^-^sf^ (3) 

In order for the solution to be regular at the origin /(O) = a(0) = 0. If the solution is to 
have finite energy, /(r) and a(r) must tend to 1 as r ^ oo. The resulting string is the well 
known Nielsen-Olesen vortex ||T^. It turns out that / and a take their asymptotic values 
everywhere outside of a small region around the string. Thus \(j)\ is constant and yl^ is pure 
gauge away from the string. The size of this region is of order r]~^. 
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Consider an extension of (jl]) to include a two-component fermion, with charge 1/2. The 
extra terms in the Lagrangian will then be 



>Cfcrmions = ipia^'D^'ilj - - [igY^lj(l)'il^'' + (h. c.)] 



(4) 



where = {—I,a^), D^il) = (^d^ — ^ieAfj^^ and = ia^ij)* is the charge conjugate of 
ij). This gives the field equations 



d.-dt 



V' - mf/(r)e*"V* 



(5) 



where the expressions (pi) and (H) have been substituted for cj) and A^, and mf = gyV- 



Following previous work by Jackiw and Rossi [11| we look for solutions with no z or t 
dependence. Such solutions will have zero energy and are called zero modes. We can see 
from (H) that they can also be taken as eigenstates of a^. The angular dependence of a 
solution satisfying a'^ip = ijj can be separated out with the ansatz 



U{r)e'^^ + V*{r)e 



(6) 



I is an arbitrary integer. We will consider the case 21 = n — \ separately. The field equations 
(I) imply 



r 2r 



dr 



n-l + l na(r)\ 
r 2r 







(7) 



Analytic solutions to these equations cannot generally be found. However it is possible to 
determine the number of physical solutions by considering their asymptotic form. 

At large r we can approximate / and a by 1. Then (0) can be solved with modified Bessel 
functions, which are asymptotically equal to e^'^f/^/r. We are interested in normalizable 
solutions, with J \^p\'^ cPx finite, so only the decaying solution is acceptable. 

When r is small, we can neglect / and a. To first order the solutions of (^) are determined 
by the angular dependence of @ 



U 
V 



, V 



U = 0{r 



.n—l\ 



(8) 



In order to match up with the one acceptable large-r solution, both the small-r solutions 
must be regular. This will only be true if < / < n — 1. This suggests there are n solutions 
of (0) if n > since there are n choices of I. In fact we have counted each solution twice, 
since putting / ^ n — 1 — / in (|6|) gives an equivalent ansatz. For every real solution of (^) 
there is also an imaginary one. This gives a total of n real solutions. 
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The above analysis needs modification when I = n — 1 — I. We can set V* = U in (^), 
giving the single equation 

I na(r)\ ^ , , 

dr-- + U + mtfU* = 9 

r 2r J 

This can be solved analytically. The solution which is well-behaved at large r is 

^(r, e) = Qr' exp J^' m,f{s) + n^ds^ (10) 

This is regular at r = if / = (n — l)/2 > 0. 

Similar analysis can be applied to the solutions satisfying a^ip = —ijj. We find a total 
of |n| solutions, which are all eigenstates of a^. Their eigenvalues are +1 if n > 0, and — 1 
if n < 0. All the solutions decay exponentially outside the string, and so are confined to it. 
They can be regarded as fermions trapped on the string. 

We can see from (^) that the solutions can easily be extended to include z and t de- 
pendence. This is achieved by multiplying ip by a{z,t), which satisfies {dz =F dt)a = 0, 
depending on whether a^ip = ±ip. Thus the trapped fermions move at the speed of light 
in the itz direction. 

These currents are conserved, and the string acts as a perfect conductor. While zero 
modes have very little cosmological significance, the lightlike currents can have dramatic 
consequences. 

3 An SO{10) GUT with Strings 

One example of a phenomenologically credible grand unified theory (GUT) has the sym- 
metry breaking 

30(10) 5^ SU{5) X Z2 

^ SU{3)cX SU{2)l xU{1)y X Z2 

^ SU{3), X U{1)q X Z2 (11) 

The discrete Z2 symmetry allows the formation of topologically stable cosmic strings. One 
possibility is an Abelian string, similar to that described in Section ||, 

*126 = e'"V„/(r) (12) 

X, = (13) 

VIO er 



00 is the usual VEV of $126, and X is a GUT gauge field broken by $126- The only stable 
Abelian strings have |n| = 1, but higher winding number strings may have long lifetimes. 
Non-Abelian strings also form in this model, but they do not have zero modes at high 
temperatures, and I will not consider them here ||n 
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The string gauge field has a nontrivial effect on the electroweak Higgs field $10 • The 
components of $10 have charges ±1/5 with respect to the generator of the GUT gauge field 
X. If n is high enough, it is energetically favourable for $10 to wind like a string. 

$10 = Hy^'K[r) + Hde-'"''hd{r) (14) 

where and Hd are the usual VEVs of the components of $io- Whether $10 winds or not, 
a nonzero electroweak gauge field is required to give a vanishing covariant derivative outside 
the string. The solution satisfies the boundary conditions /i„(oo) = hd{oo) = 6(00) = 1, 
6(0) = and /iu(0) = /id(0) = if m 7^ 0. If m = /i^ and are small, of order 



l-f^?i,cz|/l'?^ol |10l- m \s determined by the GUT string, and is equal to the nearest integer 
to n/5, so the electroweak Higgs field does not wind around a topologically stable Abelian 
string. The region of electroweak symmetry restoration is inversely proportional to the 
electroweak scale, so is far greater than the string core. 

The fermion sector of the 5*0(10) GUT contains all the Standard Model fermions, and 
an extra right-handed neutrino for each family. For simplicity I will just consider the string's 
effect on just one family, although it is easy to generalize the results. Of the three Higgs 
fields, only <I>i26 and $10 can couple to fermions. Only right-handed neutrinos couple to 

while neutrinos of either helicity couple to Hu- As the gauge symmetry unifies left- and 
right-handed particles, it is convenient to express everything in terms of left-handed spinors. 
I will use V = vl and = ic'^v'^ to express neutrino terms. Defining ^^^^ = {y'^^u)'^ , the 
resulting neutrino mass terms are 

JM I "^GUT/(r)e*"^ mu/i„(r)e^™^ \ (^) 

rriu ~ \Hu\ ~ 1 MeV is the up-quark mass and mcux ~ |</'o| ~ 10^^ GeV. Since e = 
^u/^GUT 1) the neutrino mass eigenvalues outside the string are 



yi + 4e2 + 1 _ 
iTiR — TO-GUT ;:; ~ TOqut 



2 



m-L — ttigut 7^ ~ (17j 

The mass eigenstates are then approximately + eu and v — ev'^. This illustrates the 



seesaw mechanism [13|. Although the neutrinos have the same couplings to the electroweak 
Higgs field as the up quark, the GUT Higgs ensures that is superheavy and ul is very 
light, as is required to agree with observation. Recent measurements have suggested that 
ul does indeed have a small mass . 

At high temperatures only $126 is nonzero and ttt-u = 0. The model then reduces to 
the toy model discussed in Section ^, with if) = u^. Thus |n| right-handed neutrino zero 
modes exist on the string. This implies that such strings always have conserved currents at 
high temperatures. These could allow vortons to be formed. If string loops formed at this 
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energy scale do not decay, their evolution will lead to serious conflict with observations of 
the present universe |^]. 

The situation is more complex for the neutrino fields at low temperatures since they 
couple to two Higgs fields at the same time. As in Section |2|, we will first look for solutions 
satisfying a^tp = ip. With respect to the generator of the X gauge field, ly has charge 3/10. 
It has the same charge as ^iq with respect to the Z generator. The resulting field equations 
are 

dr + -de + ly' + m^K{r)e'"'^u* + mGUT/(r)e™^i/^* = (18) 

r 2 r / 

e^' (dr+'-dg-^^+\m-'^]^-^]u + m,/ijr)e^-^- = (19) 
V r 10 r I 0} r / 

Although Jackiw and Rossi did not consider this case, it can be approached using a similar 
method to theirs. The angular dependence can be removed with the substitutions 

i.^ = (f/e*'^ + y*e^("-i-^)^) (20) 



(21) 



(The case 2Z = n — 1 will be considered later.) The resulting four complex differential 
equations are 

dr-- + U + m^K{r)W + mavTf{r)V = (22) 

n — l — l na(r)\ , , , , , , 

{r)Y + mGVTf{r)U = (23) 

m + l + l (Wm — 2n)b(r) — 3na(r)\ , , , .x 

5, + + ^ ^ W + mu/i« r t/ = 24 

r lOr / 

+ "^ + "-^ + (10^-2n)b(r)-3na(r) X ^ ^ ^ ^ 

r lOr / 

Splitting these equations into real and imaginary parts gives two identical sets of four real 
equations, and so it is only necessary to look for real solutions. 

When r is large, hu-, f, a, and b are all approximately 1. As with (|^), rearrangement of 
reduces them to modified Bessel equations at large r. The four independent solu- 
tions are proportional to 6^™^'"/-^? and 6^"^^^/-^/?, so only two of them are normalizable. 
For small r, ~ j ^ j,|n|^ g^j^^j ^j^g gauge terms are of order r. Thus we can 



ignore them when finding the leading-order terms of the small-r solutions of (22)-(|25D. 
The dominant terms of the four independent solutions are 

r , r , r , r [2b) 



If (/9 is a solution of ([2^)-(p5D for all r, which is normalisable at r = oo, then it will 
have to match some combination of the two normalizable solutions for large r. At r = 0, 
if will be made up of a combination of the solutions in (p^). So if f is to be normalizable 
everywhere, at least three of the solutions (|26|) must be well behaved at r = 0. Thus for 
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each I satisfying three of the inequahties / > 0, / < —m — 1, I < n — 1 and / > n + m, 
there will be one normalizable solution. If / satisfies all four there will be two solutions. 
Not all of these solutions are independent since the real (or imaginary) solutions for / = I' 
and I = n — 1 — I' are proportional. 

For Z = (n — l)/2 the angular dependence of (18) and (^) is removed with the substi- 
tutions 

= Ue'^^ , V = VF*e-*(™+^+')'^ (27) 
giving (after dropping gauge terms) the equations 

dr-- + U + m^K{r)W + movTf{r)U* = (28) 

r Zr / 

m + l + l (Wm — 2n)b(r) — 3na(r)\ , , . ^ 

dr + + - TTT^ ^-^)W + muhu{r)U = 29 

r lUr / 

The two real solutions have the asymptotic forms e"™'^''/ ^/r and e"*^^/ y/r, while the two 
imaginary solutions are e^^^/^/r and e~™^''/-y/r. The leading-order terms of the small-r 
solutions (real or imaginary) are 

r' , r-*"'!-' (30) 

Matching large- and small-r solutions reveals that in this case there is one real and one 
imaginary solution if < / < —m — 1. 

Combining all the above results gives a grand total of 2m (m real and m imaginary) 
normalizable solutions if m > 0, and otherwise. Surprisingly, this does not depend on 
n. A similar approach can be applied to the other components of u'^ and v to give —2m 
normalizable solutions, provided m < 0. Hence there are 2|m| possible neutrino zero modes 
after electroweak symmetry breaking. Examination of the asymptotic zero-mode solutions 
reveals that they are confined to the region of electroweak symmetry restoration. 

Since |2m| < |n|, some of the zero modes will be destroyed. For a stable n = 1 string 
all zero modes are destroyed. Thus, since higher n strings almost certainly decay, there are 
zero modes before, but not after the electroweak phase transition. The neutral current in 
the string disperses [15| and any vortons formed would dissipate after about lO"^'^ sec |^]. 
Before the electroweak phase transition from about lO^'' GeV to lO'^ GeV the universe 
would undergo a period of matter domination. Once the vortons dissipate there would be 
some reheating of the universe. However, the electroweak interactions and physics below 
the phase transition would be unaffected. 

Although we have looked at a specific 50(10) symmetry breaking, the results apply to 
most other breakings of this group, since they have the same fermion mass terms. In a 
more arbitrary C/(l)x (Standard Model) theory, the ratio of m and n could be greater than 
1, in which case extra zero modes would be created at the electroweak phase transition. 

Although we have only considered one type of model, the arguments can be applied to 
any theory with fermions coupling to cosmic strings. This has been done for a general theory 
in ref . [|l^] . The results obtained agree with simpler index theorems derived previously |l^] . 
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4 Index Theorems and Spectral Flow 

I have shown that zero modes can acquire masses at subsequent phase transitions. No 
matter how smah this mass, the spectrum of the Dirac operator changes significantly. 
If we compare the Dirac spectrum with a zero mode and a low-lying bound state with 
infinitesimal mass (Fig. |l]), we see that an arbitrarily small perturbation to the zero mode 
introduces an entire new branch to the spectrum. Any massive state gives a spectrum 
that is symmetric about both the w and k axes; there is always a reference frame in which 
the particle is at rest and others where it is moving up or down the string. Conversely, 
the zero mode, which is massless, can only move in one direction along the string and its 
spectrum is asymmetric. The transition from zero mode to low-lying bound state causes 
drastic changes in the spectrum and can be brought about by infinitesimal changes in the 
value of one Higgs field. If we consider the species with the zero mode alone, this infinite 
susceptibility to the background fields appears unphysical. However, when we include the 
massless neutrino in the 50 (10) model the spectral changes are less worrying. For a small 
coupling between the two neutrinos, both the before and after spectra have a continuum 
of massless or nearly massless states. These states can be used to build the extra branch 
of the perturbed zero-mode spectrum, allowing small changes in the overall spectrum for 
small changes in the background fields. This observation leads me to conjecture that zero 
modes can be removed only if they become mixed with other massless states. 

At the electroweak phase transition the neutrino zero mode will mix with the left-handed 
neutrino field to form a bound state. Its mass will be proportional to the neutrino mass 



inside the string, which is of order m'^^/mQ^rj, |18]. High energy currents are free to scatter 



into left-handed neutrinos off the string, so the maximal current will be very small. This 
will be insufficient to stabilize vortons. Additionally, the current will be spread over the 
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region of electroweak symmetry restoration. This is far larger than the size of a vorton, 
which is a couple of orders of magnitude greater than the GUT string radius Q. The 
current on one part of the string will then interact with current on the opposite side of the 
loop, increasing the vorton's instability. Thus vortons in SO{W) will certainly decay at the 
electroweak phase transition. 

Although I have discussed a specific type of GUT, many of the ideas I have used apply 
to a far wider range of theories. Fermion zero modes are generic in supersymmetric theo- 
ries ll^. They, too, can be destroyed at lower temperatures, in this case by supersymmetry 



breaking |2C 



5 Conclusions 

If fermions couple to a cosmic string Higgs field, their spectrum will gain extra states which 
are confined to the string core. The existence of zero-energy fermion states on strings can 
be examined analytically. Such solutions can easily be extended to give massless currents, 
which will have significant cosmological effects. In this paper I have investigated the exis- 
tence of fermion zero modes in two simple models, one of which is contained in a realistic 
5*0(10) GUT. We have seen that the microphysics of cosmic strings can be influenced by 
subsequent phase transitions. Fermion zero modes, and consequently lightlike currents on 
the strings can be created or destroyed by such phase transitions. In determining whether 
or not a cosmic string carries conserved currents it is not enough to just consider them 
at formation, but one must follow the microphysics through the multiple phase transitions 
that the system undergoes. 

It is possible for vortons formed at high energy to dissipate after a subsequent phase 
transition if the relevant fermion zero mode does not survive the phase transition, thus 
vorton bounds could be evaded. I have demonstrated this effect by analysing the neutrino 
zero modes of an SO{W) model in detail. Prior to dissipation there could be a period of 
vorton domination; after the phase transition the universe would reheat and then evolve as 
normal. 

The right-handed neutrino zero modes are removed at the electroweak phase transition 
when they mix with the Standard Model left-handed neutrinos. By considering spectral 
flow we have seen that zero modes and massless currents can only be removed by mixing 
them with another massless field. The resulting state will be a bound state or massive 
current. In the model I considered, the maximal current that the string can support after 
the electroweak phase transition is proportional to the left-handed neutrino mass. Such 
currents are far too small to stabilize vortons. 
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